Influence of initial fluctuations on geometry measures in relativistic U-U and Cu-Au 

collisions 
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In the framework of the Glauber approach we investigate the influence of the initial fluctuations 
on various measures of the initial-state geometry in 63 Cu+ 197 Au and 238 U + 238 tj relativistic ion 
collisions. Comparing variants of Glauber model (the wounded-nucleon model, the mixed model, 
and the hot-spot model) we indicate sensitivity of certain observables, in particular for the azimuthal 
eccentricity parameters as well as for the correlation of directions of the principal axes associated 
with the Fourier components. We apply GLISSANDQ in our analysis. 
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I. INTRODUCTION 

Variants of the Glauber model 0, H| , in particular the 
wounded-nucleon model 0, 0] and its extensions 0, || 
have become a basic tool in modeling the early stage of 
relativistic heavy-ion collisions. An alternative treatment 
is based on the Color Glass Condensate theory (for a re- 
cent overview see, e.g., Q and references therein). The 
Glauber model approach provides initial conditions for 
the subsequent hydrodynamic evolution. That way the 
features of the nuclei structure, such as distributions of 
nucleons in nuclei-the NN correlations Q, as well as the 
NN cross section Q show up indirectly in the measured 
observables. Moreover, as argued by Filip et al. [Tol — 
[lH, the nucleus deformation plays an important role in 
the "geometry" of the collision. Since recently collisions 
238tj + 238tj were registered at BNL RHIC, the issue of 
a proper inclusion of the nuclear deformation is impor- 
tant. Similarly, the measured collisions of asymmetric 
nuclei, 63 Cu+ 197 Au, provide valuable information of the 
initial state geometry of such systems which lead 
the geometric triangular flow [lil - flq . 

Many elements enter the modeling of the relativistic 
heavy-ion collisions: nuclear structure, models of the 
early phase, hydrodynamics, hadronization, rescattering 
of hadrons after freeze-out. Each of them brings in cer- 
tain assumptions affecting the predictions for the ob- 
served quantities. As one is primarily interested in prop- 
erties of matter created in the collision and inferred from 
the hydrodynamic phase, the early "geometric" phase 
should be modeled as accurately as possible to limit the 
uncertainty in the later stages. In particular, the effects 
of the deformation should be incorporated for such sys- 
tems as the 238 U + 238 TJ collisions, since they may lead to 



* Maciej.Rybczynski@ujk.edu.pl 
t Wojciech.Broniowski@ifj.edu, 



t Grzegorz.Stefanek@ujk.edu.pl 



r 



observable effects in the elliptic flow in most central col- 
lisions |17]. Asymmetric collisions, such as 63 Cu+ 197 Au, 
differ qualitatively from the symmetric case, since they 
give rise to the odd collective flow components not only 
from fluctuations, but also from the averaged original 
geometry of the collision, which in this case contains the 
odd Fourier components in the azimuth. 

Many recent analyses arc devoted to the studies of 
the initial geometry and its influence for the further 
stages of the evolution and the observed quantities, 
both in approaches invol ving the Glauber model [Tcj - 
l36j and the glasma [HI, l37H44j . Studies of correla- 
tions between the harmonic flow components were pre- 
sented in [H, [2^, [H, [35j]. We note that the exper- 
iment shows a strong correlation between $2 and $4 
(iH Hil, while a weak correlation between $2 and $3 
is found [47| , where $„ denotes the principal axis associ- 
ated with the n-th harmonic flow The early phase 
of the 238 U + 238 FJ collisions has been investigated in 
the Glauber approach [48M5ll ] as well as in the in AMPT 
model [Hi]. 

Notably, the generalization of the approach in the form 
of the mixed model [5|, [53|, [5J] , amending the wounded 
nucleon model with some binary collisions, leads to suc- 
cessful description of centrality dependence of multiplici- 
ties in collisions at energy range from RHIC to LHC. This 
means that we have comparable average entropy produc- 
tion from "soft" wounded nucleons and from "semi-hard" 
binary collisions. In this model the production is propor- 
tional to (1 — a)Nw/2 + aiVbin, with JVw denoting the 
number of wounded nucleons, A^in the number of binary 
collisions, and the parameter a controlling their relative 
weight. 

However, it should be stressed that there is no "unique 
Glauber model" of the early phase. The original distri- 
bution of sources (wounded nucleons, binary collisions) 
should be overlaid with some properly chosen distribu- 
tion of the entropy or energy production occurring at 
each individual NN collision. The idea goes back to the 
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TABLE I. The parameters of the Woods-Saxon potential and 
deformation coefficients for the nuclei used in our analysis, 
taken from [561 ]. 



nucleus R [fm] 


a [fm] 




ft 


63 Cu 4.206 


0.5977 


0.162 


-0.006 


197 Au 6.430 


0.45 


-0.13 


-0.03 


238 U 6.810 


0.54 


0.28 


0.093 



very beginnings of the wounded nucleon model Q , where 
a wounded nucleon produces particles with a given distri- 
bution of the average multiplicity such that the measured 
hadron multiplicity is reproduced. Adjusting the disper- 
sion of this overlaid distribution one might also reproduce 
the multiplicity fluctuations in the nucleus-nucleus colli- 
sions. Depending on the form of the overlaid distribution, 
large fluctuations of the spatial entropy distribution may 
be induced. They lead to large fluctuations of the ini- 
tial geometry, much larger compared to the naive case 
where no distribution is overlaid. These effects are stud- 
ied in detail in this work in the context of the considered 
collisions. 

A physically motivated overlaid distribution with large 
fluctuations is present in the hot-spot model, constructed 
in the spirit of [55}. The model assumes that the cross- 
section for a semi-hard binary collision producing a hot- 
spot is small, around <7bin — 2 mb, much smaller than the 
cross-section for the wounding, cr w ~ 40 — 60 mb. How- 
ever, when this rare collision occurs, it produces a very 
large amount of entropy, enhanced by the factor <7 w /<7bin 
compared to the production from a wounded nucleon. 
This factor is chosen such way that the fraction of the 
production from the binary collisions is equal (1 — a) by 
construction. 

On top of the wounded nucleons and the hot-spots one 
may still overlay a suitable statistical distribution of the 
strength of the entropy production. Here, following [l8j . 
we use the T distribution, which additionally increases 
fluctuations. The presence of hot-spots and the overlaid 
distribution, while from construction innocuous for the 
quantities dependent on the averaged densities, affects 
significantly the fluctuations, and consequently, the fluc- 
tuation measures, but also the harmonic flow measures. 
The paper is precisely focused on this class of effects. 
We study, with the help of GLISSAND0 @, the follow- 
ing aspects of the recently measured Cu+Au and U+U 
reactions at RHIC: 

• Sensitivity of the flow coefficients e„ to the in- 
creased initial fluctuations of the entropy produc- 
tion (the hot-spot + T model), as well as depen- 
dence on the nuclear deformation. Interestingly, 
we find that the fluctuations from hot-spots com- 
pletely wash-out the knee structure for the ultra- 



TABLE II. Centrality classes (ranges in RDS) for the Cu+Au 
collisions 



Centrality [%] Wounded Mixed Hot-spot 






- 5 


> 


101 


> 157.4 


> 157.7 


5 - 


- 10 


101 


- 87 


157.4 - 128.3 


157.7 - 128.2 


10 


- 20 


86 


- 61 


128.3 - 84.0 


128.2 - 84.0 


20 


- 30 


60 


- 42 


84.0 - 53.7 


84.0 - 53.5 


30 


- 40 


41 


- 28 


53.7 - 33.2 


53.5 - 32.9 


40 


- 50 


27 


- 17 


33.2 - 19.6 


32.9 - 19.7 


50 


- 60 


16 


- 10 


19.6 - 11.3 


19.7 - 11.3 


60 


- 70 


9 


- 6 


11.3 - 6.2 


11.3 - 5.9 


70 


- 80 


5 


- 4 


6.2 - 3.4 


5.9 - 3.0 



TABLE III. Centrality classes (ranges in RDS) for the U+U 
collisions 



Centrality [%] Wounded Mixed Hot-spot 






- 5 


> 


196 


> 331.3 


> 331.9 


5 ■ 


- 10 


196 


- 168 


331.3 - 272.4 


331.9 - 272.3 


10 


- 20 


167 


- 120 


272.4 - 182.9 


272.3 - 182.8 


20 


- 30 


119 


- 84 


182.9 - 119.3 


182.8 - 119.3 


30 


- 40 


83 


- 56 


119.3 - 75.4 


119.3 - 75.0 


40 


- 50 


55 


- 36 


75.4 - 44.8 


75.0 - 44.2 


50 


- 60 


35 


- 22 


44.8 - 24.4 


44.2 - 25.0 


60 


- 70 


21 


- 11 


24.4 - 12.6 


25.0 - 12.5 


70 


- 80 


10 


- 6 


12.6 - 5.5 


12.5 - 5.7 



central U+U collisions advocated in (l7j . 

• Correlations of the harmonic principal axes and 
their sensitivity to the details of the initial model 
for the Cu+Au and U+U reactions, with the con- 
clusion that the hot-spots change the qualitative 
behavior of these correlations. The effect is seen 
in the distributions of the differences of principal 
axes. 

• In our simulations we use the realistic NN collision 
profile. This profile is related to the differential 
NN cross section [57]]. We have shown previously 
that the single Gaussian wounding profile is suffi- 
ciently accurate and is the physical choice [9j], dif- 
fering in predictions for the fluctuation measures 
from the commonly used but inappropriate hard- 
sphere wounding profile. 



3 



II. MODEL 

In the analysis presented in this paper we use 
GLISSANDO 6] modified to incorporate the shape defor- 
mation of nuclei. The spatial distribution of nucleons has 
been generated according to the deformed Woods-Saxon 
density 



w 1 + exp (r - R (1 + P2Y20 + A>4o)) /a w 

The parameters used for 63 Cu, 197 Au, and 238 U nuclei 
are listed in Table |TJ In the collision, the principal axis of 
the deformed nucleus is randomly rotated in polar and 
azimuthal planes. The orientation of the deformed nu- 
cleus relative to the beam axis has a direct influence on 
centrality and eccentricity of the collisions at a given fixed 
value of impact parameter b. The effects of the deforma- 
tion are most important [Tol - fl2j for the central collisions 
of strongly deformed nuclei, such as 238 U. 

The short range NN repulsion is simulated in the 
GLISSANDO nuclear distributions via the introduction of 
nucleon-nucleon expulsion distance d. The centers of nu- 
cleons in each nucleus cannot be closer to each other than 
d, which simulates the hard-core repulsion in the nuclear 
potential. The used value d = 0.9 fm reproduces accu- 
rately [8[ the effects of the central NN correlations im- 
plemented in a more exact manner in Refs. [58l [59j. 

In the wounded- nucleon model the key entity is the 
NN collision profile, p(b), defined as the probability of 
inelastic nucleon-nucleon collision at the impact param- 
eter b. Most of the Glauber Monte Carlo generators on 
the market use, for simplicity, the hard- sphere wo unding 
profile, i.e., the collision occurs if b < ^cj w /(2tt) for the 
wounded, and b < \/<7bin/(27r) for the binary collisions. 
As shown in [57j], the NN wounding profile in the form 
of a combination of Gaussians can precisely reproduce 
the CERN ISR experimental data |60h64| on the total 
and elastic differential p + p cross section. Here we use 
a single Gaussian form, which for the studied heavy-ion 
observables is accurate enough 

p(b) = Ae-* Ab2 /°-. (2) 

The parameter A depends weekly on the collision energy 
and we use A — 0.92 for our studies @. The realistic 
Gaussian wounding profile affects important observables 
in a noticeable way, as shown in [9j. Namely, it reduces 
the eccentricity parameters as well as multiplicity fluctu- 
ations. The effects are at the level of 10-20% compared to 
the hard-sphere profile. Physically, the effects enter from 
the fact that the Gaussian profile has a tail extending to 
large values of b, thus nucleons staying far away from the 
collision center may collide with a non-zero probability. 

For various comparisons presented in this paper, we 
are using three different models as implemented in 
GLISSANDO [f|, namely the wounded- nucleon model, the 
mixed model and the hot-spot model [H, Hf| with 
overlaid T distribution. Within the Glauber approach, 



during the first stage of the collision individual inter- 
actions between the nucleons deposit transverse entropy 
(or energy). These elementary processes, stemming 
from wounded nucleons or binary collisions, are termed 
sources. A weight called relative deposited strength 
(RDS) is assigned to each source. The distribution takes 
the form 

f(RDS) = (1 - a)N w /2f w + aN hin f hia , (3) 

where a is the parameter controlling the relative weight 
of the wounded to binary sources, Nw and Abm are 
the numbers of the wounded nucleons and binary col- 
lisions, while /w and /bin are the statistical distributions 
of the RDS generated by the wounded nucleons and bi- 
nary collisions, respectively. We assume the normaliza- 
tion / duf w (u) = J duf hhl (u) = 1. 

For the wounded-nucleon model with no superimposed 
distribution a — 1 and fw(u) = 5(u — 1), i.e. the 
strengths of each source is equal. This provides a lower 
limit for the amount of initial fluctuations in the entropy- 
production mechanism. 

For the mixed model investigated here we assume a — 
0.145 (the value fitting the multiplicity distributions at 
the highest RHIC collision energy [H,[53]), and fw( u ) — 
fbin{u) = S(u — 1). This is a popular choice in many 
other simulations. 

As mentioned in the Introduction, realistically, the dis- 
tribution of sources in the transverse plane should be con- 
voluted with a statistical distribution. This convolution 
simulates the dispersion in the generated transverse en- 
tropy. Since the meaning of hot-spots is somewhat differ- 
ent in various theoretical models, let us explain in detail 
our implementation [6]. In our hot-spot model a binary 
collision is generated according to the standard criterion, 
but is accepted with the probability Obm/ow, whereas 
the RDS of a<7 w /<7bin is assigned to the hot-spot posi- 
tion in the transverse plane. We use <7bin = 2 mb, which 
means that the hot-spots occur rarely (in 5% of binary 
collisions) , but with a large weight (about 20 times larger 
than the wounded sources). The average weight of events 
is from construction equal to (1 — a)A\v/2 + aAbi n) the 
same as in the mixed model, but it can fluctuate con- 
siderably from event to event depending on how many 
hot-spots are created. In addition, we overlay the T dis- 
tribution on top of the wounded nucleons and hot-spots, 
with the distributions fw(u) = fbin(u) = r(it, k), where 

„, „ u k ~ 1 k k exp(— ku) 

r(u,«) = — '-, ue 0,oo). 4 

The r distribution gives (u) = 1 and var(it) = 1/k. In 
our simulations we use k — 2. 

We have analyzed 10 7 minimum bias 238 TJ 4- 238 u colli- 
sions and the same number of 63 Cu+ 197 Au events gener- 
ated by GLISSANDO. Since the minimum bias simulations 
contain very limited statistics of the very central colli- 
sions, we have prepared additional samples for this case. 
All events were generated with the total inelastic NN 
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FIG. 1. (Color online) The event-average harmonic eccentricities e n (n = 2, 3,4) for the Cu+Au and U+U collisions, plotted as 
function of the number of wounded nucleons, A^w- Mixed model, u w = 42 mb with the Gaussian wounding profile. The effects 
of the nuclear deformation are relevant for £2 for the most central U+U collisions. 



cross section cr w =42 mb, which is the value assumed for 
the energy y/s NN = 200 GeV. 



The RDS determines the centrality classes and its 
meaning depends on the model. The resulting central- 
ity classes are listed in Tables |TT] and IIIII From construc- 
tion, the centrality classes are essentially the same for 
the mixed and the hot-spot+r models, as the overlaid 
distributions fw{ u ) = fbin{ u ) have the mean value set 
to one. 



III. SHAPE DEFORMATION 



In this Section we investigate the Fourier components 
of the azimuthal distribution of sources in the transverse 
plane, defined as 



V(£i C cos[n(0 4 - $„)] Y + (£V r«sin[n(^ - $„)]) 2 



y ■ r n 



(5) 
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FIG. 2. (Color online) The event-average harmonic eccen- 
tricities e n (n — 2, 3, 4) for the most central U+U collisions, 
plotted as function of the relative deposited strength (RDS) 
for the mixed model and the hot-spot+F model, <r w = 42 mb 
with the Gaussian wounding profile. 



where i runs over the number of sources in each event, 
r,; is the distance of the source from the center of mass 
of the fireball, and <fii is its azimuthal angle. The weight 
associated to each source is proportional to r™ [18| , which 
makes the higher moments more sensitive to sources fur- 
ther away from the center. 

The angle $„ is adjusted in each event in such a way 
as to maximize e„, which gives a condition 



tan n<& n — 



J2i rf cos(n^) ' 



(6) 



The angle <!>„ is the phase of the principal axis of nth 
order azimuthal Fourier component. 
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centrality 20% - 50% 
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centrality 50% - 80% 




FIG. 3. (Color online) Two-dimensional event-by-event dis- 
tribution plot of £2 and £3 for the Cu+Au collisions at three 
centrality classes, the hot-spot+F model, <r w = 42 mb with 
the Gaussian wounding profile. The value of the correlation 
coefficient is given in the upper left corners. 



In Fig. [T] we test the effects of nuclear deformation. 
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FIG. 4. (Color online) Same as Fig. |3]for the U+U collisions. 



We display the dependence of event-average e n (n=2,3,4) 
on the number of wounded nucleons, N\y, in the mixed 
model. The figures show two sets of curves, solid and 
dashed, corresponding to the case with and without de- 



formation, respectively. The largest difference is ob- 
served for the eccentricity £2 for the U+U collisions at 
highest centralities. For other cases the influence of de- 
formation is small (we have checked this up to eg). 

For 63 Cu+ 197 Au we can see several rather weak effects 
connected with the deformation of the colliding nuclei. 
First, €2 is increased in most central events, decreased 
in semi-central collisions, and in addition the maximum 
value of Nw is larger. These effects can be inferred ge- 
ometrically from the fact that we collide small prolate 
63 Cu nucleus and a much bigger 197 Au nucleus with the 
oblate deformation. The triangular deformation param- 
eter £3 increases due to deformation for the mid-central 
events. 

The effects observed in the very central 
238 U + 238 U collisions of deformed nuclei are fur- 
ther investigated in Fig. [3J where we show e„ (n=2,3,4) 
as a function of RDS. We compare the mixed model and 
the hot-spot + r model in order to visualize the influence 
of the initial fluctuations. We note large qualitative 
differences between the two considered models. In the 
mixed model we see a "knee" structure, as advocated 
in [TlJ, in £2 around RDS = 400. Such a behavior is 
a consequence of the tip-tip orientation of the prolate 
238 U nuclei. In the most central events the initial 
eccentricity is reduced because the transverse profile of 
longitudinally oriented 238 U nuclei is spherical. However, 
in the hot-spot + T the knee structure disappears due 
to fluctuations from the hot-spots and the overlaid 
distribution. One thus observes that the hot-spot 
model, which assumes the production of large amount of 
entropy in rare semi- hard binary collisions (hot-spots), 
hides the purely geometrical effects in £2. The events 
from the hot-spot model exhibit much higher values 
of all e n (n=2,3,4). The increase is around 30-50% or 
even higher for the very central collisions. As shown 
in [t| [ll|, even for collisions of spherical nuclei the 
event-by-event fluctuations of eccentricity £2 are quite 
large. The deformation of colliding nuclei can slightly 
increase these fluctuations, as the similar number of 
sources can appear in collisions of nuclei with different 
orientations in the azimuthal plane. 



IV. EVENT-PLANE CORRELATIONS 

Next, we show the analysis of the correlation be- 
tween different Fourier components e„ of the azimuthal 
distribution of sources as well as the correlation be- 
tween the direction <&„ of the principal axes. Figures [3] 
and 0] present 2D correlation plots for €2 and £3 in 
three selected centrality ranges of the 63 Cu+ 197 Au and 
238 U + 238 U collisions. They show that the correlation 
between £2 and £3 is weak and does not change signifi- 
cantly with centrality. The correlation is defined in the 
standard way as 

P={(e 2 -{e 2 ))(e 2 -{e 2 ))). (7) 
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FIG. 5. The wounded-nucleon (solid) and hot-spot+F 
(dashed) model prediction for the distribution of 6($2 — $3) 
in Cu+Au collisions. 



The largest correlation (p ~ 8%) is observed in the cen- 
tral (0-20%) 63 Cu+ 197 Au events. In peripheral collisions 
(50-80%) e 2 and e 3 are slightly anti-correlated. From 
Figs. [3] and @] we note that the small correlation may be 
induced by the hnite size of the space (0 < e„ < 1) rather 
than any dynamics. 

Next, we discuss the correlation between directions $„ 
of principal axes ([6]). Most interesting is the correlation 
of the lowest rank even and odd axes, i.e., $2 and $3- A 
correlation of these axes is clearly seen in Figs. [S] and |H 
where the distribution of the angle 6($2 — $3) is shown in 



FIG. 6. Same as Fig. [5] for the U+U collisions. 

three large centrality bins for the two colliding systems 
and two models: the wounded-nucleon model and the 
hot-spot model. We note that at low centralities these 
models yield qualitatively different results. 

For the wounded-nucleon model these distributions ex- 
hibit a minimum at $2 - $3 = in central collisions (cen- 
trality 0-20%). This can be understood as follows. The 
angle $2 fluctuates around the direction perpendicular 
to the reaction plane. The triangularity angle $3 naively 
should be completely random. However, since in the col- 
liding systems there are more nucleons distributed in the 
direction parallel to the reaction plane than perpendicu- 
lar, there is a higher probability that one of the corners 
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FIG. 7. (Color online) A cartoon of the geometry of a mid- 
central collision of asymmetric nuclei. The smaller nucleus 
has a larger curvature at the geometric boundary than the 
larger nucleus, thus the overlap region has a shape of an asym- 
metric almond, hence the axes of the elliptic and triangular 
deformations form an angle of ~ 7r/6. 



of the triangle points in the direction parallel to the reac- 
tion plane. Then the relative angle between the average 
values of $ 2 and $ 3 is 7r/6 (see Fig. E}, and 6($ 2 - $3) 
shows maxima at ±7r. The effect is stronger for the asym- 
metric Cu+Au system, as in this case the difference in 
the radii of the nuclei forms on the average the arrange- 
ment of Fig. [7]). At larger centralities the nuclei collide 
more peripherally, with fewer participants. In that case 
$3 follows closer and closer $2 and the distribution of 
Figs. [5] and [5] become more sharply peaked at zero. 

In the hot-spot+r model (dashed lines in Figs. [5] and 
[5]) the situation is different. The presence of a strong hot- 
spot source collimates $2 and $3 even at low centralities, 
hence the maximum of the distribution is at $2 - $3 = 0. 

The observed qualitative difference of the distribution 
of ($2 — $3) might be used to discriminate between dif- 
ferent models of the early phase, and in particular, the 
degree of the initial fluctuations. However, we should 
keep in mind that the angles $ n are not observable. 
The observable angles are the event-plane angles, 
determined from measuring the momenta of produced 
particles. Nevertheless, it is expected that the correla- 
tions between the n = 2 and 3 axes are carried over by 
event- by- event hydrodynamics from the initial state to 
the final particle distributions 32j . The effects shown in 
Figs. [5] and |fJ] are small for the central and mid-central 
collisions (the departures from unity are at the level of 1- 
2%), but since the experimental resolution is very good, 
they should be within the experimental reach. 

A convenient measure of the above-discussed effect is 



(cos[k($ n -$ m )]), 



(8) 



where k is the least common multiple of n and m mul- 
tiplied by small natural number, and averaging is over 
the events in a given class. Such two-plane correla- 
tors are sometimes calculated for experimental data, and 
their linear combinations provide the three-plane corre- 
lators which carry additional information not accessible 
through two-plane correlators. When averaged with the 
distributions such as in Figs. [S] and |B1 the value becomes 



negative when the distributions are concave, and positive 
when they are convex. Thus the more collimated axes, 
the higher (cos[k($„ - $ m )]) is. 

The centrality dependence of (cos[k($ n — $ m )]) for 
different choices of n and m in 63 Cu+ 197 Au and 
238tj + 238tj collisions j s shown in Figs. and Fig. [9] The 
presented results suggest that the correlation between $„ 
and $ m is strong in peripheral collisions for all measured 
n and m. This simply follows from the fact that formally 
in the limit of very peripheral collisions we have just two 
sources (when there is only one source we cannot deter- 
mine the axes and the event is excluded from the sample) . 
In that case all axes are parallel to one another (and 
also e n = 1). One can see that the correlation between 
$2 and $4, as well as $ 3 and $6 is quite strong also for 
the mid-central events. The shape of the distributions 
for 63 Cu+ 197 Au and 238 U + 238 U collisions significantly 
differs only for the correlators of $ 2 and $4. We note 
a knee structure in the centrality dependence of the $2 
and $4 correlation in the most central 238 U + 238 TJ colli- 
sions for the mixed model. Similarly to the case of e 2 of 
Sec. IIIIl the structure disappears when the hot-spot+r 
model is used, showing sensitivity to the presence of the 
initial fluctuations. 

The correlation of <E> 2 and $4 in 63 Cu+ 197 Au colli- 
sions increases for central and mid-central events for 
the hot-spot model in comparison to the mixed model. 
The increase is most significant for cos(4($ 2 — $4)). 
For 238 1J -f 238 u collisions the hot-spot model does not 
change the strength of the $2 — $4 correlation but the 
knee structure observed in the most central events disap- 
pears. The correlation of $3 and $6 is increased in the 
central 63 Cu+ 197 Au and 238 U + 238 U collisions for the 
hot-spot model. 

The quantities in Figs. [8] and [9] are plotted as functions 
of the RDS. The corresponding values of centralities can 
be obtained from Tables HJ and IIIIl 



V. CONCLUSIONS 

The main point of this paper is that the initial fluctu- 
ations may qualitatively change the behavior of various 
measures of the initial geometry (eccentricity parame- 
ters, event-plane correlations). The studies of Cu+Au 
and U+U systems, recently analyzed at RHIC, clearly 
exhibit sensitivity to the choice of the selected variant of 
the Glauber-like model of the initial phase. In our study 
we have used GLISSANDD simulations with deformed nu- 
clei and the realistic Gaussian wounding profile for the 
NN collisions. We have shown that the model with a 
higher degree of fluctuations, such as the hot-spot+F 
model, may easily hide the subtle features of the initial 
geometry, such as those reflecting the deformation of the 
colliding nuclei. Conversely, a sufficiently precise mea- 
surements of the event-plane correlators may indicate a 
particular model of the early phase, provided the correla- 
tions are not largely altered by the intermediate evolution 
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FIG. 8. (Color online) The event-plane correlations measures (cos[fc($ n — <3?m)]} for the Cu+Au collisions. We compare the 
predictions of the mixed model (solid line), the wounded- nucleon model (dashed line), and the hot-spot+r model (dotted line). 
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of the system from the early phase to the freeze-out. 



fluctuations. 



In view of our results, one may reconcile the prediction 
of [ijj for the knee shape in the dependence of the eccen- 
tricity parameter on the number of produced particles in 
the most central U+U collisions. No knee in the exper- 
imental data would hint to an initial model with more 



ACKNOWLEDGMENTS 

This work was supported by the Polish Ministry of 
Science and Higher Education under Grants No. N N202 
263438, N N202 288638, and National Science Centre, 
grant DEC-2011/01/D/ST2/00772. 



[1] R. J. Glauber in Lectures in Theoretical Physics W. E. 
Brittin and L. G. Dunham eds., (Interscience, New York, 

1959) Vol. 1, p. 315 

W. Czyz and L. C. Maximon, |Annals PhysT| 52, 59 (1969) 
A. Bialas, M. Bleszyriski, and W. Czyz, Nucl. Phys. 
Bill, 46 1 (1976) 

A. Bialas, |J~T~h^s~| G35, 044053 (2008) 
D. Kharzeev and M. Nardi, |Phys. Lett.| B507, 121 (2001) 
W. Broniowski, M. Rybczynski, and P. Bozek, 
IComput. Phys. Commun~ 180, 69 (2009) 

F. Gelis, E. Iancu, J. Jalilian-Marian, and R. Venu- 
gopalan. lA~nn.Rev.Nucl.Part.Sci.l 60. 4 63 (2010) 
W. Broniowski and M. Rybczynski, Phys. Rev. C81, 
064909 (2010) 

M. Rybczynski and W. Broniowski, Phys. Rev. C84, 
064913 ( 2011) 

P. Filip, |Phys.Atom.Nucl.| 71, 1609 (2008) 

P. Filip, R. Lednicky, H. Masui, and N. Xu, |Phys.Rev.| 

C80, 054903 (2009) 

P. Filip, jNucl.Phys.P roc.Suppl.| 198, 46 (2010) 
P. Bozek, |Phys.Lett.| B7 17~287 (20 12) 

B. Alver and G. Roland, |Phys.Rev.| C81, 054905 (2010) 
B. H. Alver, C. Gombeaud, M. Luzum, and J.-Y. Olli- 
trault, |P~rys. Rev.| C82, 034913 (2010) 
H. Petersen, G.-Y. Qin, S. A. Bass, and B. Muller, 
|Phys.Rev.| C82, 041901 (2010) 
S. A. Voloshin, [Phys.Rev.Lett.l 105, 172301 (2010) 
W. Broniowski, P. Bozek, and M. Rybczynski, Phys. Rev. 
C76, 054905 (2007) 

R. Andrade, F. Grassi, Y. Hama, and W.-L. Qian, 

|J.Phys.| G37, 094043 (2010) 

P. Staig and E. Shuryak, Phys. Rev. C84. 034908 (2011) 
D. Teaney and L. Yan, |Phys.RevT~ C83, 064904 (2011) 

G. -Y. Qin , H. Petersen, S. A. Bass, and B. Muller, 
|Phys.Rev.| C82, 064903 (2010) 
J. L. Nagle and M. P. McCumber, 

044908 (2011) 

J. Xu and C. M. Ko, Phys.Rev.| C83, 
R. A. Lacey, R. Wei, N. Ajitanand, J. 
et fflZ., |Phys.Rev~| C84, 027901 (2011) 
K. Werner, I. Karpenko, T. Pierog, M. Bleicher, 
K. Mikhailov, [Phys.Rev.| C82, 044904 (2010) 
G.-Y. Qin and B. Muller, |Phy"s.Rev.| C85, 061901 (2012) 
R. S. Bhalerao, M. Luzum, and J.-Y. Ollitrault, 
|Phys.Rev.| C84, 034910 (2011) 

R. S. Bhalerao, M. Luzum, and J.-Y. Ollitrault, 
|Phys.Rev.| C84, 054901 (2011) 

F. G. Gardim, F. Grassi, M. Luzum, and J.-Y. Ollitrault, 



|Phys.Rev. C83, 

021903 (2011) 
Alexander, J. Jia, 



and 



|Phys.Rev.| C85, 024908 (2012) 
Z. Qiu, C. Shen, and U. Heinz, 



|Phys.Lett. B707, 



151 



[32 
[33 
[34 
[35 
[36 

[37 

[38 
[39 

[40 

[41 

[42 

[43 
[44 

[4& 

[46 

[47 

[48 

[49 

[50 

[51 

[52 

[53 

[54 

[55 

[56 

[57; 

[58 

[59 



(2012) 

Z. Qiu and U. W. He inz, |Phys.Rev.| C84, 024911 (2011) 
J. Xu and C. M. Ko, |Phys.Rev.| C84, 44907 (2 011) 
D. Teaney and L. Yan(201~~7 |arXiv:1206.1905 [nucl-th~ 
J. Jia and S. Mohapatra(2012), |arXiv~i~20 3.5095 [nucl-th] 
T. Hirano, P. Huovinen, K. Murase, and Y. Nara(2012), 
|arXiv: 1204.5814 [nucl-th~ 

J. L. Albacete, A. Dumitru, and Y. Nara, 
|J.Phys.ConESer~ 316, 01 2011 (2011) 
B. Muller and A. Schafer, |Piiys.Rev.|D85, 114030 (2012) 
S. Gavin and G. Moschelli, |Phys.Rev.| C85, 014905 

(2012) 

U. Heinz and J. S. Moreland, |Phys.Rev.| C84, 054905 
(2011) 

B. Schenke, P. Tribedy, and R. Venugopalan, 

|Phys.Rev.Lett~ 108, 252301 (2012) 

B. Schenke, P. Tribedy, and R. Venugopalan, Phys. Rev. 

C86, 034908 (2012) 

A. Dumitru and Y. Nara, |Phys.Rev~| C85, 034907 (2012) 
S. Gavin, G. Moschelli, and G. Moschelli, |Phys.Rev.| C86, 

034902 (2012) 

J. Adams et al. (STAR Collaboration), |Phys.Rev.Lett.| 
92, 062301 (2004) 

A. Adare et al. (PHENIX Collaboration), |Phys.Rev.Lett.| 
105, 062301 (2010) 

K. Aamodt et al. (ALICE Collaboration), |Phys.Rev. Lett" 

107, 032301 (2011) 

U. W. Heinz and A. Kuhlman, |Phys.Rev.Lett~ 94, 132301 
(2005) 

A. J. Kuhlman and U. W. Heinz, |Pl~7s~Rev~| C72, 037901 
(2005) 

H. Masui, B. Mohanty, and N. Xu, |Phys.Lett.[ B679, 440 
(2009) 

T. Hirano, P. Huovinen, and Y. Nara, |Phys.Rev-l C83, 

021902 (2011) 

M. R. Haque, Z.-W. Lin, and B. Mohanty, |Phys.Rev.[ 
C85, 034905 (2012) 

B. B. Back et al. (PHOBOS), Phys. Rev. C65, 031901 
(2002) 

B. B. Back et al. (PHOBOS), Phys. Rev. C70, 021902 
(2004) 

M. Gyulassy, D. H. Rischke, and B. Zhang, Nucl.Phys.] 
A613, 397 (1997) 

P. Moller, J. Nix, W. Myers, and W. Swiatecki, 

Atom.Data Nucl.Data Tabl. 59, 185 (1995) 

A. Bialas and A. Bzdak, Acta Phys. Polon. B38, 159 

(2007) 

M. Alvioli, H.-J. Drescher, and M. Strikman, Phys. Lett. 

B680, 225 (2009) 

M. Alvioli and M. Strikman, |Phys.Rev.| C83, 044905 



11 



I 

eg 
CO 

o 
o 



C\J 



0.6 
0.5 
0.4 
0.3 
0.2 
0.1 
0.0 



0.2 



0.1 




' 1 1 h 

u+u 

Mixed 
Wounded 
Hot-spot +r 



CO 

o 
o 

v 0.0 



CD 

i 

CM 

CD 
CO 

o 
o 



0.2 



0.1 



0.0 



0.1 



c\j 



C\J 



CO 

S o.o 




~i 1 r 



100 



200 300 
RDS 



400 



500 






(h) 



100 



200 300 
RDS 



400 500 



FIG. 9. (Color online) Same as Fig. [8] for the U+U minimum collisions. 
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